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Abstract 

We investigate the conductance and zero-frequency shot noise of interacting, multi-level quantum 
dots coupled to leads. We observe that co-tunneling assisted sequential tunneling (CAST) processes 
play a dominant role in the transition region from Coulomb blockade to sequential tunneling. 
We analyze for intermediate coupling strength the dependence of the conductance due to CAST 
processes on temperature, coupling constant, and gate voltage. Remarkably, the width of the 
CAST transport feature scales only with temperature, but not with the coupling constant. While 
the onset of inelastic co-tunneling is associated with a super-Poissonian noise, the noise is even 
stronger above the threshold for CAST processes. 

PACS numbers: 73.63.-b, 73.23.Hk, 72.70.+m 
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Introduction. Quantum dots formed in semiconductor heterostructures have been in the 
focus of interest for transport studies, because they allow considerable variation and control 
of their properties. The discrete level structure is directly observable, and the coupling to 
the electrodes can be tuned. In this letter we study a multi-level quantum dot in a regime 
of intermediate coupling strength, where besides elastic and inelastic co-tunneling the co- 
tunneling assisted sequential tunneling (CAST) processes [1, 2] play a crucial role. Transport 
signatures of CAST- processes were observed in the differential conductance [3], and their 
potential limitations for applications of devices relying on Coulomb blockade have been 
discussed. In this article we present, starting from a perturbative diagrammatic approach, 
a detailed analysis of the conductance and shot noise [4-8] in the experimentally relevant 
regime. We show that the CAST processes lead to characteristic features inside the Coulomb 
blockade region, particularly dominant in the shot noise. Different from the line width 
broadening at the threshold to sequential tunneling, the width of the CAST transport feature 
scales only with temperature, but not with the coupling constant. 

Model and technique. We consider a quantum dot with two spin-degenerate levels 
with energies = 1,2) and Coulomb interactions U and Unn for intra and inter-level 

Coulomb repulsion, respectively, 

where — cj^Cia is the number operator of electrons with spin a of the level i. The dot is 
coupled to source and drain reservoirs with a total Hamiltonian H — ifL + -f^R + -f^D + -f^T,R + 
i/T,L- The reservoirs Hr (r = L, R) are modeled by non-interacting electrons (flkor) with 
density of states pr, and the coupling to the reservoirs is described by the tunneling terms 
i^T.r with matrix elements tr (chosen independent of spin a and equal for both quantum dot 
levels). Together they define the couphng strengths = 27r|irpPr and the total hnewidth 
r = Fl -|- Fr. For transport, we apply a source-drain bias Vbias symmetrically across the 
quantum dot. In three-terminal experiments, with an additional gate voltage V^ate applied 
to the dot, the energies of the dot eigenstates are shifted by el^ateN, where N = '^i^Uia- 

We evaluate the current I = (/) and (zero-frequency) shot noise power spectrum, 
S = J^^dt{Si{t)5i{0) + 6i{0)6i{t)), with 6i{t) = i{t) - (/). The current is the 
symmetric combination, I — (7r — /l)/2 with current operators defined by — 
—i{e/h) Ylika (^r^kcrr-^^o- " h.c.^ for electrons tunneling into/out of the lead r. We proceed 
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within the frame of the real-time diagrammatic technique on the Keldysh contour [9] , which 
was extended in Ref. [5, 10, 11] to the description of the noise. The resulting current and 
shot noise are expressed in terms of transition rate matrices W (expressed in the basis of 
eigenstates of Hj)), for which analytical expressions up to second order in F^, i.e. including 
co-tunneling processes, were obtained in Ref. [12]. For fixed source-drain bias, the station- 
ary probabilities p** to occupy the quantum dot eigenstates can be determined from the 
"master" equation, Wp** = 0. Technical details on the expressions for the current and shot 
noise can be found in Ref. [5]. Rather than plotting the shot noise power spectrum S we 
will present below figures for the Fano factor, i.e. ratio between the shot noise and current, 
F = S/2eI. For uncorrelated transport events (Poissonian statistics) the Fano factor is 
F — 1. Interactions of electrons on the quantum dot may lead to bunching of electrons and 
super-Poissonian Fano factors, F > 1. In contrast, sub-Poissonian Fano factors, F < 1, 
indicate anti-bunching, which is typical for non-interacting Fermions [4] . 

Results. Transport in quantum dot systems with a large charging energy and weak 
coupling r is typically dominated by sequential processes, in which the quantum dot is 
charged by a single excess electron at a time. The applied source-drain voltage Vbias must 
provide the gate-voltagc-dcpcndcnt energy escq(V^ate) for this single-electron charging. This 
defines the threshold voltage eV^eq(V^ate) — 2eseq(V^ate) for symmetrically applied bias above 
which (at low T) sequential tunneling sets in, while states with higher excess charges are 
exponentially suppressed. This is the case inside the Coulomb blockade region, which in the 
Vgate — Hias planc takcs a diamond shape. In this regime higher order tunneling processes, 
such as co-tunneling, which transport an electron across the dot via the virtual occupation 
of dot states, dominate. Elastic co-tunneling processes, which leave the dot in its ground 
state, are the main transport mechanism at very low bias. Inelastic co-tunneling processes 
set in at a bias eV^o — ^co, where is the lowest excitation energy within the ground state 
charge sector (i.e., for fixed total charge N). Note that V^o does not depend on the gate 
voltage. 

Throughout this work we choose a temperature of ksT < 2r for which a perturbative 
expansion up to second order in the coupling constant F is appropriate [5]. We begin with 
parameters U = 20, Unn = 4, ei = —4, €2 = —2, k^T = 2r = 0.04 (all units in meV) which 
corresponds to a singly occupied ground state (N = 1) at V^ate — 0. In the left panel of Fig. 1 
we depict the differential conductance dl/dV vs. the applied bias and gate voltage. The 
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color-scale translates low conductance values with dark blue and high conductance values 
with red. The Coulomb diamond as well as two lines outside the diamond shape which 
correspond to sequential transport via excited states are clearly visible. Inside the Coulomb 
diamond, sequential transport is exponentially suppressed, while the onset of inelastic co- 
tunneling at Vco can be recognized in the figure by the horizontal step at eKo = 2 meV (for 
N = 1 we have eco — £2 — — 2meV). Note, however, that the step height depends on 
the gate voltage (see discussion below) . For a bias below that step only elastic co-tunnehng 
events contribute to the current, above both types of co-tunneling, elastic and inelastic, 
contribute. Furthermore, since for a ground state charge of N = 1 the low-energy spectrum 
satisfies the relation > 2/3escq for any gate voltage, an inelastic co-tunncling process 
occupying the upper level can directly be followed by a sequential tunneling process. Such 
transport processes are termed co-tunneling assisted sequential tunneling (CAST) processes 
and will be discussed in more detail below. Note that second order processes also play a 
significant role to the conductance at the threshold for sequential tunneling, as they lead to 
line width broadening of order F in addition to the temperature broadening already present 
for pure sequential tunneling processes [5]. 

Next, we discuss the noise properties of the system. The right panel of Fig. 1 shows 
the Fano factor in the same voltage region as the conductance plot. Outside the Coulomb 
diamond the Fano factor drops to sub-Poissonian values, 0.5 < F < 1, typical in this region. 
Inside the Coulomb diamond we observe a horizontal, gate- voltage-independent line around 
zero bias corresponding to the thermal noise in the limit ksT ^ eVuas- For <^ eVuas < 
eco the Fano factor becomes Poissonian, F = 1, since only one transport process, namely 
elastic co-tunneling, dominates. Increasing the bias further allows for inelastic co-tunneling 
beyond eVco- Now two types of co-tunneling processes with different tunneling rates compete, 
leading to a noisier current with super- Poissonian Fano factors [5]. 

From previous theoretical work [2, 12] as well as experiments [3] it is known that additional 
structure inside the Coulomb blockade regime arises if the CAST processes are only possible 
for a bias higher than the inelastic co-tunnehng energy. This is the case when the dot 
spectrum satisfies the relation < 2/3escq(V^atc)- Using parameters U = 25, f/„„ = 5, ei = 
—4, €2 = —3.5, ksT = 2r = 0.04 (all units in meV), the above relation is satisfied inside 
the Coulomb diamond with ground state charge N = 1 for the gate voltages 3.25 > eVgate > 
—0.75 as depicted in Fig. 2. Now two energy scales are visible. As before, inelastic co- 
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tunneling displays itself in form of a gate-voltage-independent conductance step (bright 
blue to dark blue) at a bias eVco — Cco = 0.5meV. The other energy scale corresponds 
to the threshold for CAST processes, eVcAST = 2[eseq(V'gate) — Cco], visuahzed by the edge 
of the bright green "angular" shape inside the Coulomb diamond. These processes are a 
combination of an inelastic co-tunneling process, exciting the quantum dot from the ground 
state to the first excited state, and a subsequent sequential tunneling processes. The angular 
edges are parallel to the Coulomb diamond edges since the gate voltage dependence of the 
CAST energy scale is the same as for the sequential tunneling threshold eVgeq — 2eseq(V^ate)- 
For both features, the value of the differential conductance depends on gate voltage. 

The angular conductance signature of CAST processes has been observed in experi- 
ments [3]. In that work transport through a multi-level quantum dot was measured and 
gate-voltage-dependent lines of the differential conductance inside the Coulomb diamond 
were frequently found. It was speculated that the visibihty of the CAST processes was due 
to (or could be enhanced by) a level-dependent asymmetry of tunneling processes. This 
would lead, however, also to an asymmetry in the gate voltage dependence, e.g. the CAST 
feature would be enhanced on the right side of the Coulomb diamond, but reduced on the 
left. When including relaxation processes of the excited states (not shown) we find that 
weak relaxation can at first enhance the CAST feature, while strong relaxation (with rate 
many times larger than F) suppresses the CAST feature such that it becomes invisible in 
the differential conductance. As has been discussed in recent work [13] such a suppression 
of CAST features could occur in many experimental situations due to relaxation processes 
involving phonons, photons or nuclear spins. However, we find that even under the influence 
of a relatively strong relaxation the Fano factor still displays the CAST processes, demon- 
strating once again the importance of shot noise measurements for a full characterization of 
the underlying transport mechanisms. The Fano factor displayed in the right panel of Fig. 2 
shows the inelastic co-tunneling excitation energy scale at eVco — 0.5 meV when it becomes 
super-Poissonian. For higher bias, the CAST angular shape inside the Coulomb diamond is 
very notable as the Fano factor reaches a yet higher value of F Ri 2, due to the additional 
competition of CAST processes with the already present inelastic co-tunneling processes. 

Analysis of CAST transport features. We now address the question how the de- 
scribed co-tunneling and CAST transport features depend on the coupling strength and 
temperature, by example of the differential conductance for the parameters of Fig. 1. The 
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left panel of Fig. 3 displays a step in the differential conductance at the inelastic co-tunnehng 
excitation energy. An analysis of the step height for different coupling strengths shows a 
linear dependence on the total coupling T. Since we have normalized the current and hence 
also the differential conductance by a factor F this means that the co-tunneling shoulder 
height scales proportional to F^, as anticipated for a second order transport process. The 
second derivative of the current depicted in the right panel of Fig. 3 further elaborates 
this fact. Here, the dotted curve for the weaker couphng F, when multiplied by the factor 
Tsoiid/rdotted — 2, practically coincides (fat dots) with the solid curve of a higher F. The 
absolute height of the co-tunneling feature is further influenced by the value of the inelastic 
co-tunneling excitation energy e^o itself, as well as the energy escq(Vgate)- There is no sim- 
ple relation for the gate- voltage dependence of the step height. This is because the second 
order transition rates depend on the energy differences between the initial/final states of 
the transition and the virtually occupied intermediate states, and there are several possible 
virtual states, namely the empty dot state and six doubly occupied states. The step height 
is lowest in the center of the Coulomb diamond, as can be seen in the left panel of Fig. 1. 

To analyze the width of the inelastic co-tunneling feature we compare the two curves 
(dashed and dotted) of the left panel of Fig. 3, where both curves correspond to the same 
total coupling F but to different temperatures. An obvious temperature broadening is 
observed. Considering the second derivative of the current for different temperatures we 
note that the width scales linearly with temperature, with full width at half maximum of 
~ 7 — 8 ksT. Though our theory can not reach the low temperature regime ksT <^ F, as 
mentioned above, the line width broadening is clearly observed at the conductance peak 
at the sequential threshold already for the value k^T = 2F used in the present analysis. 
No such broadening is observed at the inelastic co-tunneling feature, even though CAST 
processes are also involved for the considered regime [14]. This result does not support the 
interpretation of experimental flndings provided in Ref. [1]. There, a large broadening of 
a transport feature with CAST processes was interpreted as line width broadening of size 
Fr. However, the temperature dependence was not studied in that work. It remains open, 
whether there is truly a conflict of results, or merely an experimental broadening of unknown 
physical origin, which is not considered in our theory. 

Summary. In summary we discussed the influence of co-tunneling and co-tunneling 
assisted sequential tunnehng on the differential conductance and shot noise (Fano factor) 
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inside the Coulomb blockade regime. We found that the width of the features in the differ- 
ential conductance due to inelastic co-tunneling is reduced to the thermal width, even when 
CAST processes are possible, whereas the height of the feature scales with F^. Additionally, 
we have studied a situation in which the energy scale for CAST processes is separated from 
inelastic co-tunneling. The onset of CAST processes varies with gate voltage parallel to the 
Coulomb diamond edge, leading to an "angular" shape inside the Coulomb diamond. The 
Fano factor, already super-Poissonian due to inelastic co-tunneling processes, is even further 
enhanced by the CAST processes. While the conductance features due to CAST processes 
are suppressed by relaxation processes of excited states, the super-Poissonian Fano factor is 
much more robust against this decay. 

Acknowledgments. We acknowledge stimulating discussions with Maarten Wegewijs, Mar- 
tin Leijnse, Jiirgen Konig, Jonas Pedersen, and Andreas Wacker. 
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FIG. 1: (Color online) Left panel: Color-scale plot of the differential conductance (in units 10 ■ 
-ne^T/h) vs. gate voltage (V^atc) and bias voltage (Fbias)- The Coulomb diamond edges mark the 
onset of sequential tunneling. Lines outside the Coulomb diamond correspond to transport via 
excited states, whereas the horizontal step inside the diamond is due to the onset of inelastic co- 
tunneling. Right panel: Color-scale plot of the Fano factor vs. gate and bias voltage. Around zero 
bias voltage the Fano factor diverges due to thermal noise. Poissonian values, F = 1, indicate that 
elastic co-tunneling processes dominate at low bias. Super-Poissonian Fano factors occur around 
the inelastic co-tunneling excitation energy and are reduced to sub-Poissonian values outside the 
Coulomb diamond. 
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FIG. 2: (Color online) Le/t panel: Color-scale plot of the differential conductance in the region 
where Eco < — (see text). The horizontal step inside the Coulomb diamond corresponds to the 
onset of inelastic co-tunneling whereas in the green, angle shaped region additional CAST processes 
are possible. Right panel: Color-scale plot of the Fano factor vs. gate and bias voltage for the 
same parameters as on left panel. The Fano factor becomes super-Poissonian at the inelastic co- 
tunneling excitation energy scale and increases to even higher super-Poissonian values when CAST 
processes set in (edge of the black angle shaped region). 
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FIG. 3: Left panel: Differential conductance dl/dV vs. bias voltage for the same parameters as 
in Fig. 1 at fixed gate voltage Vgate = 0. The coupling F and the temperature T are varied as 
indicated on the right panel. Right panel: Second derivative of the current, d?I/dV^, vs. bias 
voltage for the same parameters. The width of the peak scales linearly with temperature. 
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